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I. INTRODUCTION 



The history of thermoelastic stresses on layered materials begins with 
Timoshenko. In a 1925 paper [Ref. 1], he first discussed the case of a directly 
bonded bi-material configuration. He studied bi-material thermostats modeled 
in the form of high aspect ratio beams subjected to uniform heating or cooling. 
Bearing stresses in closed form expressions were presented and it was noted 
that there was not an elementary way of determining the shear stress 
distribution along the bearing surface. The curvature of the beams at a given 
temperature was shown to be related to the beam geometry, the elastic 
constants of the materials, and their coefficients of thermal expansion. 

Goland and Reissner [Ref. 2] and Christensen [Ref. 3] studied the shear 
and peeling stresses along the interface of single lap adhesive joints subjected 
to axial loads. Goland and Reissner derived closed form analytical solutions 
which showed that the peak peeling and shear stresses occurred at the joint 
edges. Additionally, it was shown that all stresses depended on the joint 
dimensions and the stiffness properties of the plates. Christensen performed 
his analysis using the finite element method. He determined that a large 
stress gradient existed across the thickness of the adhesive as well as along 
the lap, especially near the ends of the lap. Neither study investigated the 
effect of temperature on system behavior. Burgreen [Ref. 4] analytically 
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computed interfacial shear stresses in layered beams and examined the 
behavior of a bi-metallic beam subjected to a transverse temperature variation. 
Lau [Ref. 5] derived closed form stress analysis solutions for a chip in a semi- 
infinite substrate and a chip on substrate with finite thickness subjected to 
temperature rises. 

A typical electronic package consists of a semi-conductor (chip) attached 
to a substrate material by an adhesive as depicted in Figure 1.1. Placing more 
power in an ever decreasing space is the trend in electronic component 
packaging. Power densities of 12.5 MW/m 3 have been reached by today's 
integrated circuit technology [Ref. 6] and future integrated circuits are 
expected to have at least twice this power density leading to higher operating 
temperatures. The resulting high operating temperatures give rise to thermal 
cycling causing thermal fatigue at the different material interfaces and can 
lead to premature failure of the circuit. To prevent failure, a maximum 
temperature of between 100°C and 110°C at a semiconductor junction has been 
set by the military [Ref. 7]. 




Figure 1.1 Electronic Package 



While conducting a literature search, numerous studies concerning the 
topic of thermal response of electronic packages was uncovered. In perhaps one 
of the most extensive studies, Suhir [Ref. 8] reviewed work conducted on bi- 
material and tri-material electronic packages that were either subjected to a 
uniform temperature or an applied axial load. For the bi-material configura- 
tion, it was determined that the maximum peeling and shear stresses along 
the interface occurred at the outer edges of the assembly. Additionally, it was 
shown that the stresses depended upon the axial and interfacial compliances 
of the component layers. 

Royce [Ref. 9], Riches [Ref. 10], and Chung [Ref. 11] all state that the 
primary cause of mechanical stresses that are induced in electronic packages 
are a result of the difference in coefficients Of thermal expansion of the 
different materials in the assembly. Royce commented that thermal cycling 
during normal operation could give rise to thermal fatigue at the various 
material interfaces and could lead to premature failure of the structure. Rich 
reported that the types of stresses that arose from thermal mismatch consisted 
of stresses normal to the chip and interfacial stresses. He predicted that 
normal stresses in the chip were at a maximum in the chip center and 
decreased to zero at the chip edge. However, the interfacial stresses were zero 
at the center and rose to a maximum at the edge. Rich concluded that 
compliance and thickness of the chip attachment materials were critical to 
controlling the stresses imposed. Two ways to reduce or avoid thermally 
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induced stresses in electronic packages was proposed by Chung. The first was 
stress reduction by coefficient of thermal expansion matching and the second 
was stress reduction by the use of extremely flexible adhesives. 

The purpose of this investigation was to gain further insight and under- 
standing into the effects that the various system parameters have on thermally 
induced behavior of an electronic package as shown on Figure 1.1. 

To achieve this, the electronic package was generalized into a tri-material 
configuration as displayed in Figure 1.2 and a parametric study was conducted 
using finite element analysis. First, the effect of the material properties, 
Young's Modulus and coefficient of thermal expansion, was examined. Next, 
the length of the midlayer was varied. Then, the effect of midlayer thickness 
on system behavior was studied. Finally, nonuniform temperature fields were 
imposed on the tri-material configuration. 




Figure 1.2 Tri-Material Configuration 
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II. FEM MODEL FORMULATION 



The purpose of this investigation was to determine the stresses resulting 
from a nonuniform temperature field on a tri-material configuration. To 
accomplish this, a finite element method (FEM) program was constructed. 
This program is presented in Appendix A and the program exec is in 
Appendix B. 

A. PROGRAM DEVELOPMENT 

A recently developed element which provides for both axial and lateral 
displacement continuity was utilized in the stress formulation. As shown in 
Figure 2. 1, each element has six degrees of freedom; axial displacements at the 
four corner points; and lateral displacements at the two ends. 

y 

A 




■o* x h 




Figure 2.1 A Typical Element With Six Degrees Of Freedom 



The lateral displacement field is determined by linear interpolation 

v t =N l Vi+N 2 vl ( 2 - 1 ) 
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where Vj and v 2 are the lateral displacements at the ends and where the linear 
shape functions Nj and N 2 are given by 

N 2 -j (2.2) 

resulting in 

2 

v(x)=J^N i (x)y j (2.3) 

i- 1 

The axial displacement field is assumed to be linear in both the axial and 
transverse direction. That is, 



(2A) 

a?=H 2 b ut + H& (2-5) 

where superscripts b and t refer to the bottom and top respectfully. The linear 
shape functions H, and H 2 are given by 

H b =H 2 b =\(h-y) H[=H 2 =l (2.6) 

h h 

resulting in 



"(jyJ-E «,(*)[#, Mu,* +H 2 (y)«/] 



(2.7) 



1*1 



From strain-displacement relations 



/ 2 2 
du di^-' b 






i*l 



( 2 . 8 ) 
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or 



and 



a (^l W i + ^2 U i ) 

»=i «r 



_du + dv 

* 




N^ut-Mph 




or 



2, f3ff. 

^5 N k“‘ 



» ««, ,1 * aw, 



The derivatives of the linear shape functions are 



i 1 



affj__ i dN 2 _ i a// 
dx l dx l dy h 



dy h 



This results in 



H i(y),..b „ K 

t ; — Oh u i) + — ; — Oh M i) 



and 



V 






V 2" V 1 



(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 

(2.14) 

(2.15) 
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Defining the axial displacement vector as 



{bJ: T =iui u{ u£ u£) (2.16) 

and the stiffness matrix as 

[K b ] = fJ^EiBFdydx (2. 17) 



where the B vector is 



, i r /dW. 

(bF=(— x -h x 

dx 1 



dN x 

~dx 



h 2 



dN 2 

~dx 







and the force vector due to temperature as 



(2.18) 



f‘f*{B}EaATdydx (2. 19) 

gives the bending matrix equation as 

( 2 . 20 ) 

Equation 2.20 defines the bending behavior. Additionally, the stiffness matrix 
can be shown as 



w-f 



2 1 - 2-1 

1 2 - 1-2 

- 2-12 1 
1-212 



( 2 . 21 ) 



Behavior due to shear is obtained as follows. Define the row vector of 
displacement degrees of freedom as 



(6/-(w, 6 u{ v, uf v 2 ) 



( 2 . 22 ) 
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The shear stiffness matrix is given by 



[K s \=f^f*{B f \G{B')cfydx 



where 



/a / dH. dHy dN, dH. 

-r± N„ 



dy 



dy dx 



dy 



[K S ]=G 



The system matrix for bending and shear is 

[*]=[jy+[*y 



dy dx 



which gives the equations for shear behavior 

[* S ]{6 S M0} 

Additionally, the stiffness matrix can be shown as 
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( 2 . 23 ) 



( 2 . 24 ) 



( 2 . 25 ) 



( 2 . 26 ) 



( 2 . 27 ) 
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Let 



Gl 


n 

IS 


Eh 


d=° 


4h 
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° 61 


2 



then 



a+2c 


-a+c d 


a -2c 


-a-c -d 


-a+c 


a+2c -d 
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a -2c d 
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-d 
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b. 



The force vector due to a temperature field is 

{F e }=Q e f o ‘f o h {B}(fydx 



where 



or 



O e =(£aA7y 



(Ft-Q'f'f* 

Jo Jo 



m 



l -H, 



dx 1 



dN, 



dx 

0 

dN, 



-H, 



— H \ 
dx 1 



dN 2 

~dx 

0 



-H r 



dydx 



(2.28) 



(2.29) 



(2.30) 



(2.31) 



(2.32) 
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Let 



then 



After solving 



A e =Q e - 

2 



(2.33) 



{F*}=A e < 



-1 

-1 

0 

1 

1 

0 



(2.34) 



[£]{6}={Fl 



(2.35) 



for {5}, the stresses can be calculated. For the e th element, the bending stress 



is 



o =Ee 

X X 



(2.36) 



or 



b e 1 , b ^4 

(«2 — 



(2.37) 



t 1 . , 5,-6, 



(2.38) 



where a x be and are the bending stresses acting on the e th element at the 
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bottom and top respectfully. The shear stress is 






(2.39) 



or 



T 



*>e 

*y 






(2.40) 



x i =G j 

where and are the shear stresses acting on the e th element at the 
bottom and top respectfully. Finally 

_gJ(62 + S 5 )-(6i + 64),2(6 6 -6 3 ) 

Xave ~ 2 ( h + l 

B. PROGRAM INPUT 

The input file for the finite element program is displayed in Appendix C. 
A mesh grid with 200 elements as shown in Figure 2.2 is generated. The input 
file allows for variation in both the element layer thicknesses and the element 
row widths. Values for Young's Modulus, Poisson’s Ratio, and the coefficient 
of thermal expansion can be inputted row by row. Additionally, the 
temperature for each level of the mesh can be inputted. Appendix C is 
annotated, showing the placement of the different inputs. All dimensions are 
in millimeters. 



(2.42). 



u [ < v 2 5 V 5 6 M . (2.41) 

h l J l h l 
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C. PROGRAM OUTPUT 



The output file from the finite element program is shown in Appendix D. 
The force and displacement vectors for each node are outputted. The top and 
bottom bending stresses and shear stress for each element are also displayed. 
Additionally, the normal stresses along the material A and material B 
interface and the material B and material C interface are outputted. 
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F igure 2.2 F mile Element Analysis Mesh F or Tri-material C onfiguration 
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III. EFFECT OF MATERIAL PROPERTIES 



This chapter studies the effect of Young's Modulus and the coefficient of 
thermal expansion on the tri-material configuration of Figure 3.1. Several 
computer runs were made assigning different Young's Modulus values to the 
three materials while maintaining a constant coefficient of thermal expansion 
value. No stresses resulted. As a consequence, the remainder of this chapter 
is concerned with variations in the coefficient of thermal expansion ( a ^). Table 
3.1 shows the parameters used for the 11 cases that were evaluated. For all 
cases, Young's Modulus was given a constant value of 100 GPa for the three 
materials and a change in temperature of 100°C was imposed. 




| ALL DIMENSIONS IN MILLIMETERS 



CENTERLINE 



Figure 3.1 Tri-material Configuration 
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TABLE 3.1 COEFFICIENT OF THERMAL EXPANSION STUDY 



CASE 


a A (lO 6 ) 


a R (lO 6 ) 


a n (10- 6 ) 


4 


100 


50 


100 


2 


100 


200 


100 


3 


300 


100 


10 


4 


300 


10 


100 


5 


100 


300 


10 


6 


150 


100 


10 


7 


200 


100 


10 


4 


400 


100 


10 


9 


200 


100 


1 


10 


200 


100 


50 


11 


200 


100 


75 






The eleven cases were grouped into four combinations to allow 
comparisons to be made. The combinations were: 

• Cases 1 and 2. a A and a c were set equal and a B was allowed to range 
from one-half a A , a c to twice a A , a c . 

• Cases 3, 4, and 5. a A , a B , and a c were given different values in different 
combinations. 

• Cases 3, 6, 7, and 8. a B and oc c were held constant but at different values 
while a A was allowed to range at values greater than a B and a c . 

• Cases 7, 9, 10, and 11. a A and a B were held constant but at different 
values while a c was allowed to range at values less than a A anda B . 
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A. CASES 1 AND 2 



Cases 1 and 2 investigated how the tri- material configuration responded 
when the top and bottom layers (materials A and C) held the same value of the 
coefficient of thermal expansion while the middle layer was given values lower 
(case 1) and higher (case 2). 

Figure 3.2 is a plot of the normal stress along the upper interface (between 
materials A and B) and lower interface (between materials B and C) as a 
function of position. For case 1, both the upper and lower interfaces are in 
tension, showing that the top and bottom layers with their larger amount of 
expansion, wanting to peel away from the middle layer. In case 2, the normal 
stresses at the upper and lower interfaces are all compressive, that is, have 
bearing stresses. The middle layer wants to expand more than the top and 
bottom layers and pushes against them. For both cases, the maximum stresses 
occur at the ends of the configuration. 

Figure 3.3 is a graph of the shear stresses acting on the upper and lower 
interfaces. The shear stresses for case 1 tend to axially extend the midlayer 
while the shear stresses for case 2 tend to axially compress the midlayer. 
Figure 3.4 shows the shear and normal stresses acting on the midlayer for 
cases 1 and 2. 



16 



100 




(edJM) ssai^s ieuu°M 



17 



Case a A (xlO /°C) otgCxlOrC) a c (xlO /°C) 

Case 1 ~ 100 50 100 

Case 2 I 100 1 200 I 100 

Figure 3.2 Normal Stress At Interlaces For Cases 1 And 2 



500 -i 1 500 




(edW) sea^greaqs 



18 



I 2 I 100 I 200 I 100 I 

Figure 3.3 Shear Stress At Interfeces For Cases 1 And 2 
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Figure 3.4 Shear And Normal Stresses On The Midlayer For Cases 1 And 2 




B. CASES 3, 4, AND 5 

Cases 3, 4, and 5 investigated the response of the tri-material configura- 
tion when all three layers had different values for the coefficients of thermal 
expansion. Case 3 was arranged with magnitudes as: a A > a B > a c , case 4 as: 
a A > a c > a B , and case 5 as: a B > a A > a c . 

Figure 3.5 is a plot of the normal stresses as a function of position along 
the upper and lower interfaces. Case 3, with the coefficients of thermal 
expansion in decreasing order from top to bottom, displayed minimal peeling 
stresses along the interfaces. Case 4, with the coefficient of thermal expansion 
lowest in the midlayer, exhibited characteristics similar to case 1. The normal 
stresses were peeling on both interfaces. Case 5, with the coefficient of 
thermal expansion in the middle layer the largest, behaved like case 2. The 
normal stresses along the upper and lower interfaces were bearing throughout. 

The plots for the shear stresses along the upper and lower interfaces are 
shown in Figure 3.6. In case 3 the shearing stresses were the lowest, with 
material B wanting to elongate axially along the upper interface. The 
midlayer wanted to axially elongate for case 4 and to axially compress for case 
5. This was the same behavior that was observed for cases 1 and 2, respectful- 
ly. Figure 3.7 displays the shear and normal stresses acting along the 
midlayer interface. 
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Figure 3.5 Normal Stress At Interfeces For Cases 3, 4, And 5 
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Figure 3.6 Shear Stress At Interfeces For Cases 3, 4, And 5 
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Figure 3.7 Shear And Normal Stresses On The Midlayer For Cases 3, 4, And 5 



C. CASES 3, 6, 7, AND 8 

This study compared the responses of the tri-material configuration when 
a B and a c were held constant at values lower than a A (with a B > a c ), and a A was 
then allowed to increase as follows: a A (case 6) < a A (case 7) < a A (case 3) < a A 
(case 8). 

Figure 3.8 is a graph of the normal stresses along the upper and lower 
interfaces for the midlayer. When the value of the coefficient of thermal 
expansion for the upper layer was kept reasonably close to the coefficient of 
thermal expansion of the midlayer (within a factor of two) as shown by cases 
6 and 7, the normal stresses tended to be minimal and were mainly bearing. 
However, when a A was further increased as in cases 3 and 8, the magnitudes 
of the normal stresses increased and became peeling throughout the upper and 
lower midlayer interfaces. 

The shear stress plot for the upper and lower interfaces is given in Figure 
3.9. Cases 6 and 7 resulted in shearing stresses at both ends of the upper and 
lower interfaces with minimal shear stresses along the majority of the 
interface lengths. In cases 3 and 8, the shearing stresses on both interfaces 
tended to axially stretch the midlayer with case 8 having a larger magnitude 
of stress. Figures 3.10 and 3.11 display the shear and normal stresses acting 
on the midlayer for cases 3, 6, 7, and 8. 
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Figure 3.8 Normal Stress At Interfaces 
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Figure 3.9 Shear Stress At Interfeces For Cases 3, 6, 7, And 8 
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a B = 100xl0' 6 /°C a c = 10xl0' 6 /°C 
E a = E b — E c = 100 GPa 

Figure 3. 10 Shear And Normal Stresses On The Midlayer For Cases 3 And 8 
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0 ^= 100xl0' b /°C a c = 10xl0' b /°C 
E A = E^=Eq= 100 GPa 

Figure 3.11 Shear And Normal Stresses On The Midlayer For Cases 6 And 7 



D. CASES 7, 9, 10, AND 11 

In this grouping, the coefficients of thermal expansion for materials A and 
B were held constant at values greater than material C (with oc A > ccg), and a c 
was allowed to increase as follows: a c (case 9) < a c (case 7) < a c (case 10) < a c 
(case 11). 

Figure 3.12 is a plot of the upper and lower interface normal stresses for 
the four cases. Similar characteristics were displayed by cases 7 and 9. Along 
the upper interface the normal stress was initially peeling, decreased to zero, 
and then became bearing. At the lower interface the normal stresses were 
initially bearing, went to zero, and then became peeling. For cases 10 and 11 
the normal stresses were peeling on both the upper and lower interfaces. 

A plot of the shear stresses acting on the midlayer is shown in Figure 
3. 13. As can be seen, along the upper interface all four cases displayed similar 
shear stress characteristics with the surface being pulled in axial tension. 
However, at the lower interface, cases 7 and 9 showed axial compression while 
cases 10 and 1 1 displayed axial tension. The shear and normal stresses acting 
on the midlayer for these four cases are shown in Figure 3.14 and 3.15. 

E. SUMMARY 

A review of the peeling and shear stresses acting along the upper and 
lower interfaces reveals force equilibrium in both the lateral and axial 
directions. However, moment equilibrium is not satisfied exactly. This 
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Figure 3.12 Normal Stresses At Interfeces For Cases 7, 9, 10, And 11 
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Figure 3.13 Shear Stress At Interfaces For Cases 7, 9, 10, And 11 
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a A = 200x 10' 6 /°C a B = 100xl0' 6 /°C 
E a = E ]3 =E c = 100 GPa 

Figure 3.14 Shear AndNormal Stresses On The Midlayer For Cases 7 And9 
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a A = 200xl0 b /°C a B = lOOxlOTC 
E a = Eg = E c — 100 GPa 

Figure 3.15 Shear And Normal Stresses On The Midlayer For Cases 10 And 11 




deviation is most probably a result of the type of finite element that was 
chosen. 

As expected, when materials A and C were given the same properties, 
symmetrical stress conditions resulted. This was clearly shown in cases 1 
and 2. When the coefficient of thermal expansion in the midlayer was smaller, 
the upper and lower layers wanted to separate from the middle layer. When 
the coefficient of thermal expansion for the midlayer was greater, the top and 
bottom layers tended to compress it. 

Maximum peeling stresses occurred when the coefficient of thermal 
expansion for the top or bottom layer was much greater than the middle layer 
as shown in case 4. Similarly, the maximum bearing stresses arose when the 
midlayer coefficient of thermal expansion was greater than the upper and 
lower layers as evidenced by case 5. 

The relative difference in the coefficient of thermal expansion between the 
layers played a similar role for the shear stresses. The maximum tendency to 
axially elongate the midlayer occurred when the coefficient of thermal 
expansion for the middle layer was much less than the upper or lower layers. 
Additionally, the maximum tendency to axially compress the midlayer 
happened when the coefficient of thermal expansion for the midlayer was much 
greater than the top or bottom layers. 
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Minimal normal stresses resulted when the values of the coefficients of 
thermal expansion for the upper, middle, and lower approached each other. 
The conditions for minimum shearing stresses were the same. 

The study of cases 7, 9, 10, and 11 produced an interesting result. As can 
be seen, cases 7 and 9 exhibited very similar stress states even though the 
coefficients of thermal expansion for the bottom layer differed by a factor of 
ten. More computer runs were conducted, further reducing the coefficient of 
expansion by several factors of ten. The changes in the state of stress was 
minimal. 

Finally, when the coefficient of thermal expansion for the tri-material 
configuration was kept constant and uniform, changing of the Young's Modulus 
values did not produce any resulting thermally induced stresses. 
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IV. EFFECT OF MIDLAYER THICKNESS 



This chapter investigates the effect of changing the thickness of the 
midlayer (material B) on the tri-material configuration. Five different 
midlayer thickness cases were studied and compared. As displayed in Figure 
4.1, the following thicknesses for material B were used: 1.00mm (case 5), 
0.50mm (case 12), 0.25mm (case 13), 0.10mm (case 14), and 0.05mm (case 15). 
For all cases, the dimensions of materials A and C were kept constant (10mm 
length and 1.0mm thickness). Additionally, the length of the midlayer was 
held at 5mm and the tri-material configuration was given a temperature 
increase of 100°C. The material properties for the five cases were given the 
following constant values: E A = E B = E c = 100 GPa, a A = 100xl0' 6 /°C, a B = 
300xl0' 6 /°C, and a c = 10xl0' 6 /°C. 

A. NORMAL STRESS DISTRIBUTION 

Plots of the normal stresses along the upper and lower interfaces of the 
midlayer for the five cases are displayed in Figure 4.2. Both the upper and 
lower normal stress distributions exhibited similar characteristics. The 
stresses were always bearing along the interfaces for all five cases with the 
maximum stresses occurring at the left and right ends of the midlayer. For 
every case, the minimum bearing stress was developed midspan of the 
interface. As the thickness of material B was decreased, the bearing stresses 



36 



Ml k 




1 mm 



~T 



Jr 



0.5 mm 






0.10 mm 





0.05 mm 



Figure 4.1 Tri-Material Configuration For Cases 5, 12, 13, 14, And 15 
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Position (mm) Position (mm) 

Case 5 Case 12 Case 13 Case 14 Case 15 

Material B Height (mm) 1 0j> 0.25 0J. 0.05 

a A = 100xl0' 6 /°C a B =300xl0' 6 /°C a c = 10xl0' 6 /°C E A = E B =E C = 100 GPa 
Figure 4.2 Normal Stresses At The Interfaces For Cases 5, 12, 13, 14, And 15 




similarly decreased. However, the shape of the normal stress distribution 
curves remained comparable. 

Comparisons were made between the maximum of the normal stresses of 
all five cases on the upper and lower interfaces and the remainder of the 
normal stresses. Figures 4.3 and 4.4 are two-dimensional and three-dimen- 
sional plots of the normal stress/maximum normal stress for the five cases 
across the upper and lower interfaces, respectfully. For both interfaces, the 
maximum normal stress occurred in case 5 at the right hand side of the 
midlayer. Similar characteristics were displayed by all cases for both the 
upper and lower. As the thickness of the midlayer decreased, the ratio of the 
normal stress to the maximum normal stress decreased. Of interest was that 
the thinner material B became, the extent of the midspan that was at a 
minimum grew. Figures 4.5 and 4.6 show the bearing stresses acting along the 
upper and lower interfaces of material B. 

B. SHEAR STRESS DISTRIBUTION 

Figure 4.7 contains plots of the shear stresses across the upper and lower 
interfaces of material B for cases 5, 12, 13, 14, and 15. Cases 5 and 12 (with 
their greater thicknesses of material B) displayed similar attributes. The 
shear stresses, upper and lower, tended to axially compress the midlayer. 
Cases 13, 14, and 15 (with midlayer thicknesses of 0.25mm, 0.10mm, and 
0.05mm) showed peak shear stresses at the left and right ends of the 
interfaces, but along the majority of the length of span the shear stresses 
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0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 



Position (mm) 

Figure 4.3 2-D And 3-D Plots Of Normal Stress/Maximum Normal Stress At The 
Upper Interface For Cases 5, 12, 13, 14, And 15 
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Figure 4.4 2-D And 3-D Plots Of Normal Stress/Maximum Normal Stress At The 
Lower Interface For Cases 5, 12, 13, 14, And 15 
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Normal Stress/Max Normal Stress 
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Case 5 Case 12 Case 13 Case 14 Case 15 

Material B Height (mm) 1 0J> 0,25 0J. 0.05 

a A = 100xl0' 6 /°C a B =300xl0' 6 /°C a 0 = 10xl0' 6 /°C E A = E B = E C = 100 GPa 
Figure 4.6 Shear And Normal Stresses On The Mkflayer For Cases 14 And 15 
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Position (mm) Position (mm) 

Case 5 Case 12 Case 13 Case 14 Case 15 

Material B Height (mm) 1 0J5 0.25 0J. 0.05 

t A = I00xl0' 6 /°C o b =300x10' 6 /°C oc Q = I0xl0’ 6 /°C E A =E B = E C = 100 GPa 
Figure 4.7 Shear Stresses At The Interfeces For Cases 5, 12, 13, 14, and 15 



fluctuated around zero. This fluctuation of positive and negative values of the 
shear stresses from point to point may have been due to a breakdown in the 
finite element program in dealing with the very small numerical values 
present. 

Comparisons of the ratio of the shear stress at a point to the maximum 
shear stress found in the group were conducted. Figures 4.8 and 4.9 are two- 
dimensional and three-dimensional plots of the shear stress/maximum shear 
stress for the five cases along the upper and lower interfaces, respectfully. For 
both the upper and lower, the maximum shear stress was located at the left 
end of material B. The two figures were similar, with the shear stress ratios 
being maximum at the extreme ends of the midlayer and then decreasing to 
a minimum towards midspan. Of interest, like the normal stress, was that as 
the thickness decreased, the range of the midspan that was at a minimum 
increased. Figures 4.5 and 4.6 display the shear stresses acting along the 
upper and lower interfaces for the midlayer. 

C. SUMMARY 

The effect of decreasing midlayer thickness on the tri-material configura- 
tion was displayed in a number of ways. The normal stresses remained 
bearing throughout as the midlayer became thinner. The maximum normal 
and shear stresses decreased with decreasing thickness. Minimal normal and 
shear stresses extended out further along the span and the stress distributions 
became more local in nature with the thinning of material B. A conclusion 
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Figure 4.8 2-D And 3-D Plots Of Shear Stress/M axnnum Shear Stress At The 
Upper Interface For Cases 5, 12, 13, 14, And 15 
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Shear Stress/Max Shear Stress 






Figure 4.9 2-D And 3-D Plots Of Shear Stress/Maximum Shear Stress At The 
Lower Interface For Cases 5, 12, 13, 14, And 15 
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Shear Stress/Max Shear Stress 





that could be drawn would be that as the thickness of the midlayer continues 
to decrease, that the only stresses of any consequence along the span exist at 
the extreme left and right ends of the interfaces. The normal stresses were 
only one-third to one-half the shear stresses. For both the upper and lower 
interfaces the maximum normal stresses remained at the far left edge in all 
cases. The maximum shear stresses along the lower interface moved from the 
left edge to the right edge as the midlayer thickness was decreased. Finally, 
the maximum shear stresses along the upper interface remained at the left 
edge for all cases. 
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V. EFFECT OF MIDLAYER LENGTH 



In this chapter, the effect of the length of the midlayer (material B) on the 
tri-material configuration was investigated. Six different midlayer lengths 
were compared and studied. As shown in Figure 5. 1, the following lengths for 
material B were used: 9mm (case 16), 7.5mm (case 17), 5mm (case 12), 2.5mm 
(case 18), 1mm (case 19), and 0.5mm (case 20). For all cases, the dimensions 
of materials A and C were kept constant (10mm length and 1mm thickness). 
Additionally, the thickness of the midlayer was held at 0.5mm and the tri- 
material configuration was given a temperature increase of 100°C. The 
material properties for the six cases were as follows: E A =E B =E C = 100 GPa, 
a A =100xl0' 6 /°C, a B = 300xl0' 6 /°C, and a c = 10xl0' 6 /°C. 

A. NORMAL STRESS DISTRIBUTION 

Figure 5.2 contains plots of the normal stresses across the upper and 
lower interfaces of material B for the six cases. To allow comparison of the 
different cases the abscissa was constructed to be a non-dimensional location, 
with the distance along the interface being divided by the total midlayer length 
for that case. 

In all six cases, the normal stresses remained bearing across the entire 
span for both the upper and lower interfaces. Cases 12, 16, 17, and 18 (with 
midlayer lengths of 5mm, 9mm, 7.5mm, and 2.5mm, respectfully) displayed 
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Case 16 L B = 9mm 




Case 17 L B = 7.5 mm 




Case 12 L B = 5mm 




Case 18 L B = 2.5mm 




Case 19 L B = lmm 




Figure 5.1 Tri-material Configurations For Cases 12, 16, 17, 18, 19, And 20 
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LOWER NORMAL STRESS (MPa) UPPER NORMAL STRESS (MPa) 





Figure 5.2 Normal Stress At Interfeces As A Function Of Non-dimensional 
Location For Cases 12, 16, 17, 18, 19, And 20 
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similar characteristics. The maximum bearing stresses occurred at the left 
and right ends of the midlayer with the minimum being developed at midspan. 
For midlayer lengths of 9mm and 7.5mm, the stress along the midspan was 
negligible. As the length was further decreased to 5mm and 2.5mm, the 
midspan normal stresses increased. The stresses at the left and right ends 
similarly increased but to a much lesser extent as the length of the above 
mentioned four cases decreased. 

Cases 19 and 20, with the short midlayer lengths of 1mm and 0.5mm, 
displayed distinctly different normal stress distributions than the four other 
cases. The left ends of the interfaces had relatively low bearing stresses. 
Across the remaining 90% of the span, the stress distribution was quite 
uniform. Figures 5.3 and 5.4 display the normal stresses acting along the 
upper and lower interfaces of the midlayer. 

B. SHEAR STRESS DISTRIBUTION 

Plots of the shear stresses along the upper and lower interfaces of the 
midlayer for cases 12, 16, 17, 18, 19, and 20 are displayed in Figure 5.5. Like 
the normal stress distribution plots, the abscissa was constructed to be a non- 
dimensional location, with the distance along the interface being divided by the 
total midlayer length for that case. 

All six cases showed similar attributes with both the upper and lower 
shear stresses tending to axially compress the midlayer. Additionally, the 
peak stresses occurred at the extreme left and right ends of the interfaces. For 
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Case 12 Case 16 Case 17 Case 18 Case 19 Case 20 

Material B Length (mm) 5 9 7J5 2J> 1 0.5 

a A = 100xl0‘ 6 /°C a B =300xl0' 6 /°C a ( ,= 10xl0 r, /°C E A = E B =E C =100 GPa 
Figure 5.3 Shear And Normal Stresses On The Midlayer For Cases 12, 16, And 17 
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Case 12 Case 16 Case 17 Case 18 Case 19 Case 20 

Material B Length (mm! 5 9 7J> 2J5 1 0.5 

a A = I00xl0' 6 /°C ccg= 300x10' 6 /°C a c = 10xl0' 6 /°C E A =E B = E C = lOOGPa 
Figure 5.4 Shear And Normal Stresses On The Midlayer For Cases 18, 19, And 20 



LOWER SHEAR STRESS (MPa) UPPER SHEAR STRESS (MPa) 





Figure 5.5 Shear Stress At Interfeces As A Function Of Non-dimensional 
Location For Cases 12, 16, 17, 18, 19, And 20 
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the greatest length of material B, case 16 with 1^= 9mm, the only non- 
negligible shear stresses appeared at the ends of the interfaces. As the 
midlayer length decreased, the extent of the span that experienced appreciable 
shear stress increased. This is clearly displayed in Figures 5.3 and 5.4. 

C. SUMMARY 

The effect of changing the midlayer length on the tri-material configura- 
tion was shown in numerous ways. As the material B length is shortened, the 
normal stresses tend to increase and become uniformly distributed across the 
span. The shear stress became distributed across a greater extent of the 
interface as the length of the midlayer lessened. The nature of the normal 
stresses was always bearing and were much less than the shear stresses 
generated. As the length of the midlayer was decreased to 1mm, shear and 
normal stresses increased. However, further reduction in length resulted in 
a slight decrease of both stresses. Except for the shortest length where it was 
at the right hand edge, the maximum shear stresses occurred at the left hand 
edge for the upper and lower interfaces. The maximum normal stresses on 
both interfaces stayed at the right edge. Finally, the maximum upper normal 
stresses were 10 to 30 percent greater than the maximum lower normal 
stresses. 
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VI. EFFECT OF NONUNIFORM TEMPERATURE DISTRIBUTION 



Chapters III, IV, and V considered the tri-material configuration as having 
a uniform temperature increase of 100°C (AT = 100°C). In this chapter the 
temperature increases of the three materials were allowed to vary, and the 
results were studied and compared. The physical dimensions of the tri- 
material configuration used in this chapter are shown in Figure 6. 1. Materials 
A and C were both given a length of 10mm and a height of 1mm. Material B 
was given a length of 2.5mm and a height of 0.1mm. The material properties 
supplied for all cases studied were as follows: E A = E B = E c = 100 GPa, a A = 
100xl0' 6 /°C, a B = 300 x 10‘ 6 /°C, and cc c = I0xl0' 6 /°C. 



10 mm 



MATERIAL A 


1 mm 




0.1 mm 

1 mm 


| : 



MATERIAL B ~ 2.5 mm ' 



Figure 6.1 Tri-material Configuration For Cases 21, 22, 23, And 24 

Four different cases were examined. In case 21, all three materials were 
allowed a temperature rise of 100°C (AT a = AT B = AT C = 100°C). A tempera- 
ture increase of 200°C for materials A, B, and C (AT a = AT B = AT C = 200°C) 
was studied in case 22. For case 23, material A was assigned a temperature 
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rise of 200°C while materials B and C were permitted 100°C temperature 
increases (AT a = 200°C, AT B = AT C = 100°C). Finally, for case 24, material A 
was allowed to have a linearly decreasing temperature profile. The top layer 
of material A was set at a temperature increase of 200°C and the temperature 
rise was then lessened until at the interface of materials A and B the 
temperature increase was 100°C. This profile is shown in Figure 6.2. 
Materials B and C were allowed 100°C rises in temperature (AT a = linear, AT B 
= AT C = 100°C). 

200 °C 



150 °C 

125 °C 
113 °C 
106 °C 
100 °C 

Figure 6.2 Material A Linear Temperature Distribution For Case 24 

A. NORMAL STRESS DISTRIBUTION 

The normal stresses acting along the upper and lower interfaces of the 
midlayer for the four cases are shown in Figure 6.3. Similar characteristics 
were displayed for both the upper and lower normal stress distributions. The 
maximum stresses occurred at the far left and right ends of the interfaces and 
were always bearing. For the majority of the midspan the normal stresses 
were minimal. 
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LOWER NORMAL STRESS (MPa) UPPER NORMAL STRESS (MPa) 
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Figure 6.3 Normal Stresses At The Interfaces For Cases 2 1, 22, 23, And 24 
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A comparison between case 21 and case 22 revealed an interesting point. 
Case 22, with a tri- material configuration temperature rise of 200°C, resulted 
in normal stresses which were twice the value of the normal stresses of case 
21 (AT a = AT b = AT C = 100°C). This showed the direct relationship between 
system temperature increase and the resulting increase in normal stresses. 

When only material A was given a larger temperature rise, as in case 23 
(AT a = 200°C, AT B = AT C = 100°C), the change in system behavior for normal 
stresses when compared to a uniform temperature field (case 2 1) was minimal. 
Along the upper interface the response was nearly identical, with the only 
difference being the marginally larger bearing stress at the extreme right end 
for case 23. Across the lower interface, case 23 normal stress values were 
continuously slightly less than those of case 21. 

When the temperature rise in material A was allowed to linearly decrease 
through its thickness, as in case 24, the resulting normal stress distributions 
were the lowest of the four cases studied. The maximum stresses developed 
were approximately one-half of the maximum for the uniform temperature field 
case. 

Figures 6.4 and 6.5 display the normal stress distributions acting on the 
upper and lower interfaces of the midlayer for the four cases. 

B. SHEAR STRESS DISTRIBUTION 

Figure 6.6 contains the plots of the shear stress distributions across the 
upper and lower interfaces for cases 21, 22, 23, and 24. The resulting curves 
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Figure 6.4 Shear And Normal Stresses On The Midlayer For Cases 21 And 22 
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Figure 6.5 Shear And Normal Stresses On The Midlayer For Cases 23 And 24 
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Figure 6.6 Shear Stresses At The Interfeces For Cases 21, 22, 23, And 24 
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were similar for all four cases. The maximum shear stresses occurred at the 
extreme right and left ends of the interfaces. Along the majority of the length 
of span the shear stresses fluctuated around zero. This fluctuation may have 
been due to a breakdown in the finite element program in dealing with the 
very small numerical values present. 

As with the normal stresses, a comparison of case 2 1 (AT a = AT B = AT C = 
100°C) with case 23 (AT a = AT B = AT C = 200°C) shows that the shear stresses 
at every point for the higher temperature increase were twice the value of 
those stresses for the lower temperature rise. This also displayed the direct 
relationship between system temperature increase and the resulting increase 
in shear stresses. 

Case 23 (AT a = 200°C, AT B = AT C = 100°C) and case 24 (AT a = linear, AT B 
= AT C = 100°C) produced shear stress distribution values very close to case 21. 
This suggests that for the particular tri-material configuration chosen for the 
chapter, increasing the temperature rise for one material does not appreciably 
change the system’s shear stress response. 

Displays of the shear stress distributions acting on the upper and lower 
interfaces for the four cases are shown in Figures 6.4 and 6.5. 

C. SUMMARY 

A study of the effects of a nonuniform temperature field produced a 
number of interesting results. A doubling of the system temperature rise 
resulted in a doubling of the normal and shear stresses. Normal stresses 
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remained bearing in different uniform and nonuniform temperature fields. 
The normal and shear stress distributions maintained the same general shapes 
regardless of temperature distribution. Increasing the temperature rise in just 
one material resulted in minimal change in system response. Finally, a 
temperature rise which was allowed to linearly decrease through the thickness 
of one material resulted in the lowest normal stress distributions. 



65 



VII. SUMMARY AND CONCLUSIONS 



The results of this investigation into the behavior of a tri-material medium 
subjected to a temperature field leads to several conclusions. These 
conclusions are broken down into the following topics: effect of material 
properties, effect of thickness of midlayer, effect of length of midlayer, and 
effect of nonuniform temperature field. 

A. EFFECT OF MATERIAL PROPERTIES 

When both the temperature rise and the coefficients of thermal expansion 
were kept uniform throughout the system and only Young's Modulus was 
allowed to vary, no normal or shear stresses were developed. 

Changing the coefficient of thermal expansion produced several results. 
As the mismatch between the coefficients of thermal expansion of the materials 
became greater, the normal and shear stress magnitudes along the upper and 
lower interfaces became greater. The shear stresses developed along the 
interfaces were much greater than the normal stresses. Depending on the 
values of the coefficients of thermal expansion, the normal stresses were either 
bearing or peeling. Finally, along the upper and lower interfaces the normal 
and shear stresses were very local in nature. Only in the immediate area of 
the left and right corners of the interface were the stresses significant. 
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B. EFFECT OF MIDLAYER THICKNESS 



For the particular tri-material configuration studied the normal stresses 
remained bearing as the midlayer became thinner. The maximum normal and 
shear stresses decreased with decreasing material B thickness. As the 
midlayer was made thinner, the normal and shear stress distributions became 
more and more local in nature. It is believed that as the midlayer becomes 
extremely thin, the only stresses of any consequence would exist at the 
extreme left and right corners of the interfaces. The shear stresses generated 
were two to three times greater than the normal stresses. As the thickness of 
material B was decreased, the location of the maximum shear stress along the 
lower interface moved from the inner edge to the outer edge. However, along 
the upper interface the maximum shear stress remained fixed at the inner 
edge. Lastly, the maximum normal stress stayed at the left edge for both the 
upper and lower interfaces. 

C. EFFECT OF MIDLAYER LENGTH 

With the tri-material configuration used the resulting shear stresses were 
much greater than the normal stresses. The normal stresses developed were 
always bearing. As midlayer length was decreased to 1mm, the normal and 
shear stresses increased. Further shortening of the length of the midlayer 
resulted in slight decrease of the shear and normal stresses. On both 
interfaces the maximum shear stress occurred at the left hand edge. The only 
exception was the shortest length, where the maximum shear stress appeared 
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at the right hand edge. The maximum normal stress was experienced at the 
right hand edge for both the upper and lower interfaces. Additionally, the 
upper maximum normal stress was approximately 10 to 30 percent larger than 
the lower maximum normal stress. 

D. EFFECT OF NONUNIFORM TEMPERATURE DISTRIBUTION 

When the temperature rise for all three materials of the system was 
doubled, the result was a doubling of the normal and shear stresses. The 
normal stresses remained bearing and the shear and normal stress distribu- 
tions maintained the same general shapes regardless of the temperature 
distribution. Increasing the temperature rise in just the upper layer resulted 
in minimal change in system response. Finally, a linearly decreasing 
temperature rise gave the lowest magnitude normal stress distributions. 

E. RECOMMENDATIONS FOR FURTHER RESEARCH 

A number of different areas for further research in this topic exist. 
Validate the tri-material response to nonuniform temperature fields by using 
a larger number of system configurations. Alter the finite element code to 
allow study of the system response to transient temperature fields. Refine the 
tri-material configuration to more closely match electronic component 
packaging. Finally, investigate the effect of thermal cycling on fatigue. 
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APPENDIX A 
FEM PROGRAM 

**** ****************************************************************** 

* 

FEM PROGRAM FOR WELD PROBLEM * 

* 

A******************************************************************** 

COMPILE AND RUN THIS CODE WITH THE SINGLE EXEC COMMAND: WELD 



NEL 

NODES 

NDOF 

NBC 

IBC 

NFC 

IFC 

I COR R 

NMAT 

IMAT 

ITYPE 

E 

G 

X 

HITE 

ZLEN 

EK 

SYSMAT 

FLOAD 

FORCE 

U 

DELTMP 

ALPHA 

NNSCOL 

NSCOL 

N1ROW 

N3ROW 

PLSTRS 

N3STRT 



THE NUMBER OF ELEMENTS 

THE NUMBER OF SYSTEM NODES 

THE NUMBER OF SYSTEM DEGREES OF FREEDOM 

THE NUMBER OF ESSENTIAL BOUNDARY CONDITIONS 

THE SUPPRESSED DEGREES OF FREEDOM 

THE NUMBER OF APPLIED FORCES 

THE SYSTEM NODES WITH APPLIED LOADS 

THE CORRESPONDENCE BETWEEN LOCAL AND SYSTEM DOF 

NUMBER OF DIFFERENT MATERIALS 

A PARTICULAR MATERIAL 

* 0 FOR B.C.'S AT LEFT END ONLY, 

= 1 FOR B.C.'S AT BOTH LEFT AND RIGHT ENDS 

YOUNGS MODULUS OF IMAT MATERIAL 

SHEAR MODULUS OF IMAT MATERIAL 

THE X COORDINATE OF A NODE 

THE HEIGHT OF AN ELEMENT 

THE LENGTH OF AN ELEMENT 

THE 6 BY 6 ELEMENT STIFNESS MATRIX 

THE NDOF BY NDOF SYSTEM STIFFNESS MATRIX 

THE MAGNITUDE OF AN APPLIED LOAD 

THE SYSTEM FORCE VECTOR 

THE SYSTEM DISPLACEMENT VECTOR 

THE TEMPERATURE CHANGE 

THERMAL COEFFICIENT OF EXPANSION 

THE NUMBER OF NONSOLDERED COLUMNS 

THE NUMBER OF SOLDERED COLUMNS 

THE NUMBER OF ROWS IN THE CHIP 

THE NUMBER OF ROWS IN THE SUBSTRATE 

THE PEELING STRESS AT THE NODE TO THE RIGHT OF THE 

COLUMN NUMBER 

ELEMENT NUMBER BEFORE START OF LOWER SOLDER 



INCLUDE 'WELDCOMN' 

CHARACTER* 20 NAME 
COMMON /WORKSP/RWKSP 
REAL RWKSP (137047) 

CALL IWKIN( 137047) 

* 

READ *, NAME 

* 

* NOTE; INPUT THE NAME OF THE FILE TO BE OPENED AS FOLLOWS: 

* '/FN FT' 

* THE FILENAME MUST BE PRECEDED WITH A '/ AND THE FILETYPE MUST BE 

* FOLLOWED BY A THE ONLY SPACE IS BETWEEN FN AND FT. 

* 

OPEN (9, FILE* NAME ) 

* 

* INPUT THE FEM PARAMETERS 

* 

READ (9,*) NEL, NODES, NMAT, NDOF, NBC, NFC, ITYPE 
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READ ( 9 , *) (IBC(I), I « 1 , NBC) 



( 



DO 50 IB * 1 , NBC 
50 CONTINUE 

IF (NFC .NE. 0) THEN 

READ (9 , *) ( IFC ( I ) , FLOAD(I) , I - 1,NFC) 

DO 70 IF - 1 , NFC 
70 CONTINUE 
ELSE 

END IF 

DO 80 IM - 1 , NMAT 

READ (9 , 81 ) E ( IM) , POISS(IM), ALPHA (IM) , HITE (IM) 

80 CONTINUE 

81 FORMAT (4F10.7) 



READ (9 , 22 ) (X(IN), IN - 1, NODES) 

DO 999 ILL - 1, NEL 

READ (9 , 20) IL,IMAT(ILL) , ( ICORR ( ILL, IN) , IN - 1,6) , INODE ( ILL, 1 ) , 

& INODE (ILL, 2) 

999 CONTINUE 
20 FORMAT ( 1014 ) 

22 FORMAT (7F10. 5) 

READ (9,*) DELTMP 

READ (9 , *) (TEMP (I) , I - 1,NMAT+1) 

WRITE ( 14 , *) ' ' 

WRITE (14,*) 9 DELTMP * ', DELTMP 

* 

* CALCULATE THE SHEAR MODULI I 

* 

WRITE (14,*) ' ' 

WRITE ( 14 , *) 9 IMAT# , E, G MODULI I, ALFA, POISSON AND HEIGHT FOLLOW 9 

WRITE ( 14 , * ) 9 9 

DO 100 IM « 1, NMAT 

G ( IM) - 5 . *E ( IM) / ( 12 . * ( 1+POISS ( IM) ) ) 

WRITE (14, 105) IM, E(IM), G(IM), ALPHA (IM) , POISS(IM), HITE(IM) 

100 CONTINUE 

105 FORMAT (3X, 13 , 2X, 5E12 . 5) 

* 

WRITE ( 14 , *) 9 9 



DO 350 I * 1 , NDOF 
FORCE (I) ■ 0.0 
DO 340 J * 1 , NDOF 
SYSMAT ( I , J) * 0.0 
340 CONTINUE 
350 CONTINUE 

* 

* THE ELEMENT DO LOOP TO 500 FORMS THE SYSTEM MATRIX AND VECTOR 

* 

DO 500 IEL « 1 ,NEL 
IM - IMAT(IEL) 

* 

* FORM THE ELEMENT STIFFNESS MATRICES 

* 

IL - INODE (IEL, 1) 

IR * INODE (IEL, 2 ) 
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ZLEN(IEL) * X ( IR) - X(IL) 
H - HITE(IH) 

EE - E(IM) 

GG * G ( IM) 

A * GG*ZLEN (I EL) / (4 . *H) 

B - GG*H/ZLEN (IEL) 

C « EE*H/ ( 6 . *ZLEN ( IEL) ) 

D «= 0 . 5*GG 



AP2C - 


A ♦ 2 . *C 


ZMAPC - 


- 


A + C 


AM2C « 


A - 


- 2.*C 


ZMAMC - 


- 


A - C 


EK ( 1 , 1 ) 


b 


AP2C 


EK ( 1 , 2 ) 


S 


ZMAPC 


EK ( 1 , 3 ) 


m 


D 


EK (1,4) 


b 


AM2C 


EK ( 1 , 5 ) 


B 


ZMAMC 


EK ( 1 , 6 ) 


m 


-D 


EK ( 2 , 1 ) 


m 


ZMAPC 


EK ( 2 , 2 ) 


B 


AP2C 


EK(2,3) 


S 


-D 


EK (2 , 4) 


S 


ZMAMC 


EK ( 2 , 5 ) 


S 


AM2C 


EK ( 2 , 6 ) 


B 


D 


EK ( 3 , 1 ) 


B 


D 


EK ( 3 , 2 ) 


B 


-D 


EK ( 3 , 3 ) 


B 


B 


EK ( 3 , 4 ) 


B 


D 


EK ( 3 , 5 ) 


B 


-D 


EK ( 3 , 6 ) 


B 


-B 


EK ( 4 , 1 ) 


B 


AM2C 


EK ( 4 , 2 ) 


S 


ZMAMC 


EK ( 4 , 3 ) 


S 


D 


EK (4,4) 


= 


AP2C 


EK ( 4 , 5 ) 


= 


ZMAPC 


EK ( 4 , 6 ) 


= 


-D 


EK ( 5 , 1 ) 


B 


ZMAMC 


EK ( 5 , 2 ) 


* 


AM2t 


EK ( 5 , 3 ) 


B 


-D 


EK (5,4) 


3= 


ZMAPC 


EK ( 5 , 5 ) 


S 


AP2C 


EK ( 5 , 6 ) 


* 


D 


EK ( 6 , 1 ) 


S 


-D 


EK ( 6 , 2 ) 


B 


D 


EK ( 6 , 3 ) 


= 


-B 


EK ( 6 , 4 ) 


= 


-D 


EK ( 6 , 5 ) 


*= 


D 


EK ( 6 , 6) 


S 


B 



* ASSEMBLE ELEMENT MATRICES INTO GLOBAL STIFFNES MATRIX 

* 

DO 480 ID * 1,6 
IN * ICORR( IEL, ID) 

DO 480 JD * 1,6 
JN « ICORR( IEL, JD) 

SYSMAT(IN, JN) * SYSMAT ( IN , JN) + EK(ID,JD) 

480 CONTINUE 
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< 500 CONTINUE 



CONSTRUCT THE FORCE VECTOR DUE TO TEMPERATURE CHANGE 

WRITE ( 14 , * ) ' ' 

IF (DELTMP .NE. 0.0) THEN 
DO 550 IEL - 1, NEL 
IM « IMAT (IEL) 

CHI2 « E( IM) * ALPHA ( IM) * (2 . *TEMP (IM) +TEMP ( IM+1) ) *HITE ( IM) / 6 . 
CHI 1 - E ( IM) * ALPHA ( IM) * (TEMP ( IM) +2 . *TEMP ( IM+1) ) *HITE ( IM) / 6 . 

11- ICORR(IEL,l) 

12- ICORR (IEL, 2 ) 

14- ICORR( IEL, 4 ) 

15- ICORR (I EL, 5) 



550 



FORCE (II) 
FORCE (12) 
FORCE (14) 
FORCE (15) 
CONTINUE 
ELSE 
ENDIF 



FORCE (II) 
FORCE (12) 
FORCE (14) 
FORCE (15) 



CHI1 
CHI 2 
CHI1 
CHI 2 



MODIFY STIFFNESS MATRIX AND FORCE VECTOR FOR BOUNDARY CONDITIONS 

DO 610 I - 1 , NBC 
II - IBC(I) 

DO 600 J » 1 , NDOF 
SYSMAT (II , J) « 0.0 
600 CONTINUE 

SYSMAT ( II , II ) - 1. 

FORCE ( II ) - 0.0 
610 CONTINUE 

IF (ITYPE .EQ. 1) THEN 
READ (9,*) ILASTTOP 
WRITE ( 14 , * ) ' ' 

WRITE (14,*) 'B.C. S AT BOTH ENDS; ILASTTOP = 

DO 620 IN - 1 , NDOF 
SYSMAT (ILASTTOP, IN) - 0.0 
CONTINUE 

SYSMAT ( ILASTTOP , ILASTTOP) - 1. 

SYSMAT (ILASTTOP, NDOF) « -1. 

FORCE (ILASTTOP) - 0.0 
ELSE 
ENDIF 



620 



ILASTTOP 



630 



IF (NFC .NE. 0) THEN 
DO 630 IF - 1 , NFC 
II - IFC(IF) 

FORCE (II) - FORCE (II) + FLOAD(IF) 
CONTINUE 
ELSE 
ENDIF 



SOLVE THE SYSTEM OF EQUATIONS FOR THE DISPLACEMENTS 
CALL DLSARG (NDOF, SYSMAT, NPARAM, FORCE, 1, U) 
PRINT THE FORCE AND DISPLACEMENT VECTORS 
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WRITE (14,*) ' ' 

, WRITE (14,*) ' FORCE AND DISPLACEMENT VECTORS FOLLOW' 

WRITE( 14 , *) ' ' 

DO 680 I DOF « 1 , NDOF 

WRITE (14,641) IDOF , FORCE(IDOF), U(IDOF) 

680 CONTINUE 

641 FORMAT (2X, 14, 4X, E12.5, 10X, E12.5) 

* 

* CALCULATE THE STRESSES 



DO 700 IEL « 1, NEL 

11 « IC0RR(IEL, 1) 

12 « ICORR (IEL, 2 ) 

13 « ICORR (IEL, 3 ) 

14 - ICORR (IEL, 4 ) 

15 « ICORR (IEL, 5) 

16 « ICORR (IEL, 6) 

U1BOT * U(I1) 

U1TOP « U ( 12 ) 

VI « U(I3) 

U2BOT « U ( 14 ) 

U2TOP - U (15) 

V2 - U(I6) 

IM - I MAT (IEL) 

EE - E ( IM) 

GG « G ( IM) 

SIGMABOT ( IEL) * EE*((U2BOT - U1BOT) / ZLEN ( IEL) - ALPHA (IM) * 

& TEMP ( IM+1) ) 

SIGMATOP (IEL) * EE*((U2TOP - U1TOP) /ZLEN ( IEL) - ALPHA (IM) *TEMP ( IM) ) 
TAUELEM (IEL) « GG*( ( (U1TOP -U1BOT) + (U2TOP - U2BOT) ) / 

& ( 2 . *HITE (IM) ) + (V2 - VI) /ZLEN (IEL) ) 

700 CONTINUE 



* 



WRITE ( 14 , * ) 
WRITE ( 14 , *) 
WRITE ( 14 , * ) 
WRITE ( 14 , * ) 



t t 

' STRESSES FOLLOW' 
/ / 



SGMABQT 



SGMATOP TAUELEM ' 



PSI « E(1)*ALPHA(1)*DELTMP 
DO 800 IEL - 1 , NEL 
IM « I MAT (IEL) 

SIGBTNDL - SIGMABOT (IEL) /PSI 
SIGTPNDL « SIGMATOP (IEL) /PSI 
TAUNDL « TAUELEM (I EL) /PSI 

WRITE (14,801) IEL, SIGMABOT ( IEL) , SIGMATOP (IEL) , TAUELEM (IEL) 



800 CONTINUE 

801 FORMAT ( (2X,I4,6(4X,G12.5) ) ) 
WRITE ( 14 , *) ' ' 



READ (9,*) NNSCOL, NSCOL, N1ROW, N3ROW 

* 



NCOL « NNSCOL + NSCOL 

* 



DO 1300 IEL « 6, 70 
IM * I MAT (IEL) 

PLFCE(IEL) « TAUELEM (I EL) *HITE(IM) 
1300 CONTINUE 
II - 6 
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DO 1500 ICOL - 1,8 
* I ELF « 6 + (ICOL - 1) 

I ELL * I ELF + 64 
PLF - 0.0 

DO 1400 I EL - I ELF , I ELL, 14 

PLF - PLF + PLFCE ( IEL+1) - PLFCE(IEL) 

1400 CONTINUE 

UPLSTR ( ICOL) - 2. ♦PLF/ (X(II+2) - X (II ) ) 

WRITE (14,31) ICOL, UPLSTR(ICOL) 

II - II + 1 
1500 CONTINUE 

31 FORMAT ( IX , 12 , ' UPPER PEEL STRESS - ',G12.5) 

* 

WRITE (14, * ) ' ' 

DO 1600 I EL - 121, 197 
IM - IMAT(IEL) 

PLFCE (IEL) - TAUELEM ( IEL) ♦HITE ( IM) 

1600 CONTINUE 
II - 6 

DO 1800 ICOL - 1,8 
IELF - 121 + (ICOL - 1) 

I ELL - IELF + 76 
PLF - 0.0 

DO 1700 IEL - IELF, IELL, 17 

PLF - PLF + PLFCE ( IEL+1) - PLFCE (IEL) 

1700 CONTINUE 

BPLSTR (ICOL) «(2.^PLF/(X(II+2) - X(II) ) ) ♦ (-1. 0) 

WRITE (14,32) ICOL, BPLSTR ( ICOL) 

II - II + 1 
1800 CONTINUE 

32 FORMAT (IX, 12, ' LOWER PEEL STRESS - / ,G12.5) 

* 

♦ COMPARE MOMENTS DUE TO BENDING AND PEEL STRESSES 

* 

AVEFCE ( 1 ) « 0 . 5* ( SIGMATOP ( 1) + SIGMABOT ( 1) ) ♦HITE ( 1) 
ZMOMARM ( 1) - .0562 - .014 

AVEFCE ( 2 ) - 0 . 5 ♦ ( SIGMATOP (15) + SIGMABOT ( 15) ) ♦HITE (2 ) 
ZMOMARM ( 2 ) - .0562 - .028 - .007 

AVEFCE ( 3 ) - 0.5* (SIGMATOP (29) + SIGMABOT(29) ) ♦HITE(3) 
ZMOMARM ( 3 ) * .0562 - .028 - .014 - .0035 

AVEFCE (4 ) - 0 . 5 ♦ ( SIGMATOP (43) + SIGMABOT (4 3 )) ♦HITE ( 4 ) 
ZMOMARM ( 4 ) - .0035 + .00175 

AVEFCE (5) - 0.5* (SIGMATOP (57) + SIGMABOT (57) ) ♦HITE(4) 
ZMOMARM (5) - .00175 
ZMOMBDG - 0.0 
DO 1900 1-1,5 

ZMOMBDG - ZMOMBDG + ZMOMARM(I) ♦AVEFCE (I) 

1900 CONTINUE 

* 

ZMOMPEEL - 0.0 
DO 2000 I - 1,8 
II - 7 + I 
JJ - II - 2 
KK « II - 1 

W - 0 • 5* (X ( II ) - X(JJ)) 

Y - 0.25*(X(JJ) + X(II) +2 . ♦X (KK) ) 

ZMOMPEEL - ZMOMPEEL + W* Y^UPLSTR ( I ) 

2000 CONTINUE 

C WRITE ( 14, ♦) 'MOMPEEL - ', ZMOMPEEL 
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< ZMOMBDG « 0.0 
IM « 0 
Y2 - 0.056 
DO 3000 I - 1,57,14 
IM - IMAT(I) 

Y1 « Y2 - HITE (IM) 

AOV3 - ( SIGMATOP ( I ) -SIGMABOT (I) ) / ( 3 . *HITE ( IM) ) 

BOV2 « 0 . 5* ( SIGMATOP (I) - (Y2/HITE(IM) ) * (SIGMATOP ( I ) -SIGMA BOT ( I ) ) ) 
ZMOMBDG - ZMOMBDG + AOV3* ( Y2**3-Y1**3) +BOV2* ( Y2**2-Y1**2 ) 

Y2 - Y 1 
3000 CONTINUE 

SIGMAT - 1.052 
SIGMAB - 2.3012 
Y2 - .001 
Y 1 - .0005 

B = SIGMAT - (Y2/ . 0005) * ( SIGMAT- SIGMAB) 

A = ( SIGMAT-SIGMAB) / ( . 0005) 

STOP 

END 
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APPENDIX B 
PROGRAM EXEC 

FORTVS2 WELD (AD(DBL4) ICA) 

FILEDEF 14 DISK WELD OUTPUT (RECFM FB LRECL 102 BLOCK 10200) 
LOAD WELD (START 
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APPENDIX C 
INPUT FILE 



1 < 



""100. 


E+09 




.3 


100 


.E-' 


06 


• 


500 


100. 


E+09 






.3 


100 


.E- 


06 


• 
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E+09 






.3 
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.E-' 


06 
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.E- 


06 
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100. 


E+09 
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.E- 


06 


.0625 
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.3 
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.E- 


06 


. 


010 
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.3 
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.E~ 


06 
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2 < 


.3 
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.E- 


06 


>3 • 


040 
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. E- 
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.E- 


06 


• 
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.3 


10 


.E- 
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• 
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2. 
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6. Layer Temperature Rise 
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APPENDIX D 
OUTPUT FILE 

DELTMP - 200.000000000000000 



I MAT# , E, G MODULI I , ALFA, POISSON AND HEIGHT FOLLOW 



1 


0.10000E+12 


0 . 32051E+11 


0. 10000E-03 


0. 30000E+00 


0. 50000E+00 


2 


0.10000E+12 


0 • 32051E+11 


0.10000E-03 


0. 30000E+00 


0. 25000E+00 


3 


0 • 10000E+12 


0 . 3205 1E+11 


0. 10000E-03 


0. 30000E+00 


0 . 12 500E+00 


4 


0 • 10000E+12 


0.32051E+11 


0.10000E-03 


0.30000E+00 


0 . 62 500E-01 


5 


0 . 10000E+12 


0.32051E+11 


0. 10000E-03 


0 . 30000E+00 


0.62500E-01 


6 


0.10000E+12 


0 . 32 051E+11 


0 . 3 0000E-03 


0 . 30000E+00 


0.10000E-01 


7 


0 . 10000E+12 


0 . 3205 1E+11 


0.30000E-03 


0 . 30000E+00 


0.20000E-01 


8 


0 . 10000E+12 


0 . 32051E+11 


0 . 30000E-03 


0 . 30000E+00 


0.40000E-01 


9 


0. 10000E+12 


0.32051E+11 


0.30000E-03 


0 . 30000E+00 


0 . 2 0000E-01 


10 


0 . 10000E+12 


0.3205 1E+11 


0 . 30000E-03 


0 . 30000E+00 


0. 10000E-01 


11 


0 . 10000E+12 


0 . 32051E+11 


0. 10000E-04 


0. 30000E+00 


0 . 62 500E-01 


12 


0 . 10000E+12 


0. 32051E+11 


0. 10000E-04 


0 . 3 0000E+00 


0. 62500E-01 


13 


0 • 10000E+12 


0 . 32051E+11 


0 . 10000E-04 


0. 30000E+00 


0.12500E+00 


14 


0 . 10000E+12 


0 . 32 051E+11 


0.10000E-04 


0 . 30000E+00 


0 . 2 5000E+00 


15 


0 • 10000E+12 


0 • 3 2051E+1 1 


0.10000E-04 


0. 30000E+00 


0 . 50000E+00 



FORCE 


AND DISPLACEMENT VECTORS 


FOLLOW 


1 


0 . OOOOOE+OO 


0. 00000E+00 


2 


0 . 00000E+00 


0. 00000E+00 


3 


0 • OOOOOE+OO 


0. 00000E+00 


4 


0. OOOOOE+OO 


0. 00000E+00 


5 


0. 00000E+00 


0. OOOOOE+OO 
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6 


0 . 00000E+00 


0.00000E+00 


7 


0 . 00000E+00 


0 . OOOOOE+OO 


8 


0 . 00000E+00 


0.16022E-01 


9 


0 . 00000E+00 


-0.51672E-02 


10 


0 . 00000E+00 


0 . 12042E-01 


11 


0 . 00000E+00 


0.10033E-01 


12 


0 . 00000E+00 


0.90521E-02 


13 


0 . 00000E+00 


0 . 85658E-02 


14 


0 . 00000E+00 


0 . 803 3 IE-02 


15 


0 . 00000E+00 


0 • 3 2451E-01 


16 


0 . 00000E+00 


-0 . 2 1458E-01 


17 


0 . 00000E+00 


0 . 24 079E-01 


18 


0 . 00000E+00 


0 - 19888E-01 


19 


0 . 00000E+00 


0 • 17 771E-01 


20 


0 • 00000E+00 


0.167 16E-01 


21 


0 • 00000E+00 


0.15765E-01 


22 


0 • 00000E+00 


0 • 49296E-01 


23 


0 . 00000E+00 


-0 . 50058E-01 


24 


0 . 00000E+00 


0 . 3 6110E-01 


25 


0 . 00000E+00 


0.29507E-01 


26 


0 . 00000E+00 


0 • 26255E-01 


27 


0 • 00000E+00 


0 . 24626E-01 


28 


0 . 00000E+00 


0.22845E-01 


29 


0 . 00000E+00 


0 . 57798E-01 


30 


0 . 00000E+00 


-0.69213E-01 


31 


0 . 00000E+00 


0 . 42 198E-01 


32 


0 . 00000E+00 


0.34279E-01 


33 


0 . 00000E+00 


0 . 30238E-01 


34 


0 . 00000E+00 


0 . 28255E-01 


35 


0 . 00000E+00 


0 . 26546E-01 
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36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 



0 . OOOOOE+OO 


0.61934E-01 


0 . OOOOOE+OO 


-0 . 8018 IE-01 


0. 00000E+00 


0.45342E-01 


0. OOOOOE+OO 


0.36811E-01 


0. OOOOOE+OO 


0. 32342 E- 01 


0. OOOOOE+OO 


0 . 2994 3E-01 


0 . 15000E+08 


0.27144E-01 


0. OOOOOE+OO 


0. 65005E-01 


0. OOOOOE+OO 


-0.88399E-01 


0. OOOOOE+OO 


0 • 47763E-01 


0. OOOOOE+OO 


0.38867E-01 


0. OOOOOE+OO 


0.34269E-01 


0. OOOOOE+OO 


0 . 31904E-01 


0. OOOOOE+OO 


0.29455E-01 


0. OOOOOE+OO 


0 . 68161E-01 


0. OOOOOE+OO 


-0.97019E-01 


0. OOOOOE+OO 


0.50215E-01 


0. OOOOOE+OO 


0.41041E-01 


0. OOOOOE+OO 


0 . 36432E-01 


0. OOOOOE+OO 


0.34170E-01 


0. OOOOOE+OO 


0. 32225 E- 01 


0. OOOOOE+OO 


0.74620E-01 


0. OOOOOE+OO 


-0 . 11546E+00 


0. OOOOOE+OO 


0 . 55219E-01 


0. OOOOOE+OO 


0.45459E-01 


0. OOOOOE+OO 


0. 40623E-01 


0. OOOOOE+OO 


0.38212E-01 


0. OOOOOE+OO 


0.35508E-01 


0. OOOOOE+OO 


0 . 87846E-01 


0. OOOOOE+OO 


-0 . 157 15E+00 
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66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 



0. 00000E+00 


0 . 65382E-0 1 


0 . 00000E+00 


0 . 54 158E-01 


0 . 00000E+00 


0.48522E-01 


0. 00000E+00 


0.457 04E-01 


0. 00000E+00 


0.43179E-01 


0 . 00000E+00 


0. 10123E+00 


0 . 00000E+00 


-0.20527E+00 


0 . 00000E+00 


0.75574E-01 


0 * 00000E+00 


0. 62744E-01 


0 . 00000E+00 


0.56355E-01 


0 . 00000E+00 


0 . 53 160E-01 


0 . 00000E+00 


0 . 49681E-01 


0 . 00000E+00 


0.11433E+00 


0 • 00000E+00 


-0 . 2 5976E+00 


0 . 00000E+00 


0.85778E-01 


0 . 00000E+00 


0.71538E-01 


0 . 00000E+00 


0.64345E-01 


0 . 00000E+00 


0.60737E-01 


0. 00000E+00 


0 . 57389E-01 


0 . 00000E+00 


0 . 12051E+00 


0 . 00000E+00 


-0.28934E+00 


0 . 00000E+00 


0.90866E-01 


0 . 00000E+00 


0.76222E-01 


0, 00000E+00 


0.68790E-01 


0 . 00000E+00 


0.64988E-01 


0 . 00000E+00 


0.60910E-01 


0 . 00000E+00 


0 . 12 339E+00 


0 . 00000E+00 


-0.30467E+00 


0 . 00000E+00 


0.93387E-01 


0 . 00000E+00 


0.78682E-01 
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96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 

108 
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-0 . 13 044E+09 
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0 . 44885E+07 
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-0. 16488E+06 
-0. 11367E+08 
-0 . 44518E+08 
-0 . 84443E+08 
-0. 12074E+09 
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0. 14019E+08 
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0 . 81349E+08 
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